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Dissociation of superoxide dismutase 1 dimers is enhanced by glutathionylation, although the dissociation 
constants reported to date are imprecise. We have quantified the discreet dissociation constants for 
wild-type superoxide dismutase 1 and six naturally occurring sequence variants, in their unmodified and 
glutathionylated forms, at the ratios expressed. Unmodified superoxide dismutase 1 variants that shared 
similar dissociation constants with SODl WT had disproportionately increased dissociation constants when 
glutathionylated. This defines a key role for glutathionylation in superoxide dismutase 1 associated familial 
amyotrophic lateral sclerosis. 

Approximately one in five cases of familial amyotrophic lateral sclerosis (fALS) is due to naturally occurring 
mutations in the gene encoding superoxide dismutase 1 (SOD1) 1 . Over 150 such mutations have been 
reported to occur, and these mutations can occur at virtually any residue of the primary sequence. Many of 
these mutations confer a significant level of instability upon the native dimer, which may ultimately unfold and 
form insoluble cytoplasmic aggregates - a hallmark of SOD1 associated fALS pathology 2 . Recently, evidence has 
emerged to suggest that SODl WT aggregation is also involved in some cases of the considerably more prevalent 
sporadic ALS 3 " 5 . 

The prevailing hypothesis of SOD1 aggregation suggests that dissociation of the native homodimer is a key 
step, as SOD1 monomers have been identified as a common misfolding intermediate 6 ' 7 . This mechanism is 
supported by the extraordinary difference in dissociation constant (K D ) reported for the SODl A4V dimer 
(3 uM) 8 , compared with SODl WT (2.25 nM) 9 ; i.e. the dissociation of SODl A4V into monomers is approximately 
1000 times more facile at concentrations found in vivo. Furthermore, chemical cross -linking of the aggregation 
prone variant, SODl A4V can prevent its aggregation in vitro*, while cross-linking of other mutants conferred 
greater thermostability 10 . 

A compounding factor in SOD1 dissociation is the post-translational modification (PTM) of Cyslll by the 
ubiquitous antioxidant, glutathione (GSH), which has been demonstrated as being the dominant PTM in human 
tissue 9 (Fig. la; Fig. Sla). Similar to many residues in the primary chain of SOD1, Cyslll 11 is proximal to the 
dimer interface and, not surprisingly, the introduction of GSH has been proposed to augment dissociation 11 . We 
explored this hypothesis for a set of human SOD1 variants with well-defined and broad ranging rates of disease 
progression; from 1.2 years survival post-onset for SODl A4V , to approximately 17.6 years for SODl H46R . 

Results 

Our initial assay, using mass spectrometry, revealed that all variants were glutathionylated to some extent when 
produced in E. coli, from a low of 3% for SODl vl48G to 36% for SODl A4V (Fig. la). When these data were arranged 
in order of disease progression, a plausible correlation was evident for SODl A4V , SODl G93A , SODl G37R and 
SODl H46R ; however substantial glutathionylation of non-pathogenic SODl WT compared to the aggressive variant 
SODl vl48G , suggested that this modification alone is not predictive of disease progression. Indeed, while the 
sequence of these SOD1 variants may determine their susceptibility to glutathionylation, further work is required 
to confirm this in ALS patients. 
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To compare the relative dissociation propensity of our SOD1 var- 
iants, we performed nanoelectrospray-ionization mass spectrometry, 
with conditions optimized to preserve solution phase structures dur- 
ing ionization and subsequent in vacuo detection 12 (Fig. lb; Fig. S2). 
Each of the SOD1 variant spectra exhibited peaks arising from two 
dominant ion populations at —2900 m/z and —3200 m/z, corres- 
ponding to the 10 + and 11 + ions of their respective dimers. Only 
minor peaks were observed for the 9 + and 12 + charge states, indi- 
cating that the dimers had negligible variation in their tertiary and 
quaternary structures. On closer inspection, a series of three maxima 
were identified within the 10 + and 11 + peaks, arising from dimers 
composed of; (i) two unmodified SOD1 (t/SODl) monomers; (ii) 
wSODl paired with gsSODl; (iii) two gsSODl monomers. With 
the exception of SODl A4V , monomers were significantly less abund- 
ant than dimers, but similarly confined to two (7 + and 8 + ) charge 
states (Fig. lb & inset). For all variants, the most abundant monomer 
observed wasgsSODl (Fig. lc), even though wSODl constituted the 
majority of total protein. 

The mixture of SOD1 dimers was further interrogated using tan- 
dem mass spectrometry to isolate an ion population exclusive to the 
dimer (D 11+ ), prior to collision-induced dissociation (CID). At the 
minimum accelerating voltage (6 V) required to traverse an argon - 
filled CID region, the spectrum contained a single large peak of 
isolated D 11+ ions (Fig. S3a). At an accelerating potential of 20 V, 
two additional peaks arose at —2700 m/z and 3200 m/z, being due to 
the dissociation of D 11+ into charge-complementary M 5+ and M 6+ 
monomers. These peaks increased in intensity at 40 V and 60 V, 
after which the more asymmetric dissociation pathway via M 4+ 
and M 7+ monomer ions became competitive. This hierarchical gas- 
phase molecular division, based on charge, was also observed for 
^sSODl (Fig. S3b). We compared the M 6+ abundance for wSODl 
andgsSODl under low (10 V) versus high (70 V) energy CID con- 
ditions: at 10 V, £sSODl A4V , #sSODl G93A and gsSODl G37R had a 
greater proportion of M 6+ than wSODl indicating that less energy 
is required to remove gsSODl monomer from the dimer than 
wSODl. A comparison of the spectra at 70 V versus 10 V shows that 
gsSODliwSODl is larger at low energy CID for all seven proteins, 
which supports the hypothesis that modification by glutathione 
destabilizes the dimer. 

A caveat to the gas -phase analysis of molecular interactions is the 
absence of bulk water which, to varying degrees, can affect the 
strength of non-covalent association between proteins 13 . Thus, it 
was possible that the more facile removal ofgsSODl from the gas- 
phase dimer (Fig. lb&c) is not necessarily reflective of the aqueous 
environment found in vivo. To address this uncertainty, we prepared 
dilution series of the SOD1 variants ranging from 20 uM down to 
625 nM, prior to mass analysis using identical instrument conditions 
for all samples. The spectra acquired at each concentration provided 
a snapshot of all species present in solution and so, by plotting the 
proportion of monomer versus concentration (Fig. 2a; Fig. S4), we 
were able to calculate discreet dissociation constants for (wSODl) 2 
and (gsSODl) 2 (Table 1; Fig. 2b). Unlike the homodimers, which 
dissociate into their respective and unique monomer pairs, products 
of the mixed dimer (wSODligsSODl) contributed to both popula- 
tions of monomer, and were not so readily estimated. In an effort to 
account for monomers originating from the mixed dimer, we mon- 
itored the decline in (wSODligsSODl) 11+ ion intensity with success- 
ive dilution, and subtracted an equivalent quantity from each 
monomer signal. Although imprecise, we believe this correction atte- 
nuated, to some extent, the confounding effects of the mixed dimer 
species. 

All wSODl variants had K D values greater than wSODl WT (9 ± 
1 nM), implying that these single residue substitutions are suf- 
ficient to disturb the quaternary structure (Table 1). The most 
substantial augmentation was observed for wSODl A4V (1.1 ± 
0.1 uM), being approximately 100-fold greater than wSODl WT , 
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Figure 1 | Mass spectrometry analysis of SOD1 variants, (a) Percentage of 
recombinantly expressed SOD1 that was glutathionylated (yellow), as 
determined by the area under the peaks of MS spectra, compared to disease 
progression (grey), expressed in truncated life years as log e [(75 — (25 + 
d)] y where d equals the reported average disease duration, assumed life 
expectancy is 75 years and age at onset of disease is set to 50 years, (b) Mass 
spectrum of native SODl WT (main) and SODl A4V (faint) with species 
charge states labeled (D = dimer, M = monomer). Additional maxima 
within the peaks are the result of glutathionylation of one or both subunits 
to give a mixed (*) or gsSODl (**) dimer. Inset shows the ratio of dimer to 
monomer for the two proteins, derived from charge state peak areas in the 
spectrum, (c) Expansion of the M 7+ region showing peaks corresponding 
to uSODl (red) andgsSODl (yellow), (d) The percentage of gsSODl 
monomer (yellow) compared to disease progression (grey) for all SOD1 
variants. The gsSODl measurements were performed in triplicate and 
error bars represent standard error of the mean. 
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Figure 2 | SOD1 monomer abundance increases with dilution, (a) of SOD1 variants showing changes in the proportion of monomeric uSODl (red) and 
gsSODl (blue) as a function of protein concentration: (i) SOD WT , (ii) SODl A4V , (iii) SODl G37R , (iv) SODl H46R , (v) SOD1 D90A , (vi) SODl G93A , (vii) 
SODl vl48G . Analyses were performed in triplicate at each concentration and error bars represent standard deviation about the mean, (b) Relationship 
between disease duration and K D for uSODl (red) and gsSODl (blue). 



whereas uSODl G93A (12 ± 1 nM) had a K D only marginally higher 
than wSODl WT . Significantly, the K D value of each gsSODl species 
was higher than the corresponding uSODl, with the largest (13- 
fold) increase observed for #sSODl G93A , followed by#sSODl A4V (9- 
fold),£sSODl G37R (6-fold) and£sSODl WT (4-fold) (Table 1). There 
was insufficient #sSODl vl48G and £sSODl D90A to accurately cal- 
culate their K D values. 

Discussion 

Glutathionylation of SOD1 expressed in E. coli, has previously been 
reported to occur at each of the four cysteine residues (Cys6, Cys57, 
Cyslll, and Cysl46) as detected by mass spectrometry 14 . Contra- 
stingly, our analysis of recombinant SODl WT demonstrated that no 
more than one glutathione was bound per subunit, and this was the 
case for all SOD1 variants in this study (Fig. SI). This discrepancy 



stems from the protocol used by Furukawa and coworkers 14 , which 
directed expression of SOD1 into inclusion bodies or soluble apo- 
proteins followed by refolding and metal enrichment after extraction 
from bacterial cells 15 . In the present work, SOD1 was co-expressed 
with the yeast-copper-chaperone for SOD1 (yCCS) in the presence of 
copper and zinc, with subsequent purification of soluble SOD1 from 
the bacterial cytosol 16 . The yCCS facilitates folding of the nascent 
chain, and formation of the disulfide bond between Cys57 and 
Cysl46 14 . This disulfide effectively locks the P -barrel structure, ren- 
dering these Cys residues inaccessible to larger molecules; ipso facto, 
unlikely sites of glutathionylation. Furthermore, it has been demon- 
strated that Cys6 is buried inside the P-barrel of catalytically active 
SOD1, i.e., inaccessible to aqueous reductants 17 . We assessed the 
dismutase activity of our variants using in-gel zymography, which 
yielded positive results for all but the metal -binding site mutant, 



Table 1 SOD1 dissociation constants as determined by native mass spectrometry. The dissociation constants are means ± 
from 3 distinct data sets with each containing 5 data points (see supplementary figure 4) 
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SOD1 H46R ( Fig Slb ) Therefore, although we did not explicitly dem- 
onstrate Cysl 1 1 to be the site of our observed glutathionylation, it is 
inferable that the only cysteine residue sufficiently exposed and 
therefore substantially modified by glutathione was Cysl 11. 

It is well documented that patients suffering from SOD 1 -assoc- 
iated fALS have variable survival times, dependent upon the particu- 
lar amino acid mutation 18 . Considering the role played by dimer 
dissociation, as an entree to aggregation, we scrutinized our data 
for evidence of a link between patient survival and iC D (Fig. 2b). 
For wSODl, the K D values were poorly associated with disease dura- 
tion (R 2 = 0.58), primarily due to the outlier characteristics of 
wSODl G93A . In contrast, iC D values forgsSODl showed a much stron- 
ger correlation (R 2 = 0.79), indicating that glutathionylation may be 
responsible for the disparity between perceived SOD1 stability and 
observed patient survival. Consider the G93A mutation. This is 
reported as having little effect upon dimer dissociation 19 , and our 
results for the K D of wSOD 1 G93A confirm this (Fig. 3a); yet the survival 
time of patients carrying this mutation is, on average, a mere 2.2 
years. Our ability to distinguish between, and calculate separate K D 
values for wSODl andgsSODl resolves this apparent inconsistency - 
the iC D of#sSODl G93A is 13 times higher than that of wSODl (Fig. 3b; 
Table 1) - highlighting the significance of this methodology. 

Our data can be used to reframe the original hypothesis in more 
specific terms: i.e., it is not only (the extent of) SOD1 glutathionyla- 
tion that influences the rate of disease progression, but also the iC D of 
the glutathionylated species. Glutathionylation, however, is not a 
requirement for pathogenicity, as sequence modification alone of 



SOD1 vi48G causes a 3o_fbld increase in K D over SODl WT . Similarly, 
our extracted K D value for wSOD 1 A4V reflects the considerable instab- 
ility of this non-glutathionylated species. These mutations are 
located deep in the inter-subunit cleft, where small molecular volume 
changes are able to disrupt the dimer interface leading to more facile 
dissociation 20 (Fig. S5). 

In contrast to Vall48 and Ala4, Gly37 and Gly93 are spatially 
remote to the dimer interface, being situated in the P3/P4 Greek 
key loop and p5/p6 inter-strand loop respectively (Fig. S5). It is 
notable that these loops are proximal to each other and the Leu38 
cork, while Cysl 1 1 is within the P6/P7 Greek key loop at the opposite 
end of the P -barrel. Prior to glutathionylation, the individual muta- 
tions G37R and G93A impart negligible effect upon the K D of SOD1; 
i.e., in isolation, these mutations are not pathologically destabilizing 
(Table 1). Rather, it is the synergy of glutathionylation and point 
mutation at the P-barrel poles' that accounts for the significantly 
higher K D values reported in this study, and presumably, the malig- 
nance of these variants in vivo. It stands to reason then, that glu- 
tathionylation of SOD1 alters the K D of a multitude of individual 
variants in a unique, allosteric-like fashion, and that it is the com- 
bination of K D (MSOD1) and K D ^ sSOD1) (Fig. 3c) which impact the patho- 
genicity of the protein. Therefore, we propose an elaboration to the 
contemporary notion in which SOD1 dissociation is the primary 
underlying cause of SOD1 fALS. 

Specifically, dissociation of the SOD1 dimer is a biophysical con- 
sequence of either, (i) a potent destabilizing mutation such as A4V, 
or (ii) glutathionylation of a weak or non- destabilizing mutation 
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Figure 3 | The profound effect of glutathionylation on SOD1 dissociation constants. Histograms comparing K D values of (a) wSODl, (b) gsSODl, 
and (c) uSODl + gsSODl, to disease progression (grey), (d) A scheme of the proposed factors affecting SOD1 dimer dissociation. Briefly, dimer 
dissociation occurs due to mutation, glutathionylation, and cellular stresses. The increased monomer population is more likely to unfold and increase 
cellular stress. Glutathione (GSH) levels in the cell are elevated during periods of oxidative stress, thereby increasing the chances of SOD1 becoming 
glutathionylated. 
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such as G93A. Consequently, there is a twofold increase in the 
abundance of monomers, which are more susceptible to misfold into 
toxic conformations compared to the dimer 6,7 (Fig. 3d). If the cellular 
clearance rate of misfolded SOD1 is inadequate, increasingly fre- 
quent hydrophobic interactions ultimately give way to aggregation, 
forming larger amorphous structures that stress the cellular 
machinery 21 . Many neurodegenerative diseases, including ALS are 
also associated with oxidative stress 22 , which in turn can lead to an 
up-regulation of glutathione 23 . Indeed, an increase in GSH has been 
observed in ALS patient CSF 24 . Analysis of human tissue from ALS 
patients, and in cell experiments, are currently under way to further 
investigate the role of glutathionylation in SOD1 associated fALS. 

Methods 

Materials. All solutions were prepared using Milli-Q ultra-purified water using the 
Millipore purification system (Massachusetts, USA). All materials were of analytical 
grade. Tryptone, Tris-base, DTT, f3-mercaptoethanol, bromophenol blue, brilliant 
blue, ammonium persulphate, SDS, ampicillin sodium salt, and IPTG were from 
Amresco (Ohio, USA). NaCl, ammonium sulphate, agar, acetonitrile, Tris-HCl, 
tetramethylethylenediamine, ammonium acetate and RNase were from Sigma- 
Aldrich (Missouri, USA). DNase was from Roche Diagnostics (NSW, Australia). 
Carbenicillin was from Carbenicillin Direct (UK). Copper sulphate pentahydrate was 
from Ace Chemical Company (SA, Australia). Zinc sulphate heptahydrate was from 
Hopkin and Williams LTD (Birmingham, UK). Formic acid and acetic acid were from 
Univar (NSW, Australia). Methanol was from Ajax Finechem (NSW, Australia). Gold 
coated borosilicate electrospray capillaries were made in-house. 

Expression and purification of recombinant SOD1. Plasmids encoding SODl WT , 
SODl A4V , SODl G37R , SODl H46R , SOD1 D90A , SODl G93A , and SODl vl48G were a gift from 
the Oliveberg Group (Sweden). Protein expression and purification were performed 
according to 16 . Briefly, the SOD1 plasmids were transformed into chemically 
competent BL2 1 * E. coli using heat shock. Following transformation, the protein was 
expressed in the presence of copper and zinc ions with the addition of IPTG. The 
expressed SOD1 was purified via size exclusion chromatography (Hiload 16/60 
Superdex 75 PG, GE USA) and anion exchange chromatography (Hiscreen Capto-Q, 
GE USA). 

Sample preparation for mass spectrometry. SOD1 samples were desalted and buffer 
exchanged into 200 mM NH 4 OAc (pH 7), using gel filtration chromatography 
(Superdex 200 10/300 GL, GE USA). Briefly, SOD1 samples were concentrated to 
20 mg/ml and 100 ul was loaded onto the column at a flowrate of 0.4 mg/ml. The 
concentration of SOD1 was determined at 265 nm using a UV/Vis 
spectrophotometer, and a molar extinction coefficient of 18 700 M _1 cm _1 25 . 

Mass spectrometry. SOD1 samples in 200 mM NH 4 OAc were denatured using 
formic acid and acetonitrile with the final concentration of each being 10% and 40% 
respectively. Mass analyses were performed on a Q-ToF ULTIMA mass spectrometer 
(Waters, UK) in positive ion mode with nanoelectrospray source. Instrument 
parameters included: capillary 1.5 kV, sample cone 200 V, desolvation gas flow 
180 L/h. Native protein analyses were performed in 200 mM NH 4 OAc on a SYNAPT 
G2 HDMS (Waters, UK) in positive ion mode. Instrument parameters were similar to 
Ruotolo et al. (2008) 26 with the exceptions: capillary voltage 1.5 kV, sample cone 
30 V, cone gas 70 L/h, trap collision voltage 6 V, ion-transfer stage pressure 4.15 e" 1 
mbar, ToF analyzer pressure 2.38 e" 6 mbar. All SOD1 variants were analyzed at a 
concentration of 20 uM. For MS/MS experiments, the trap collision energy was 
increased from 10 to 100 V in increments of 10 V with 30 scans at each increment. 
The dilution series mass spectra were acquired at concentrations of 20, 15, 10, 7.5, 5, 
2.5, 1.25 and 0.625 uM with the following instrument parameters: capillary voltage 
1.52 kV, sample cone 137 V, cone gas 78 L/h, trap collision voltage 16 V, transfer CE 
8 V, ion transfer stage pressure 3.72 e" 4 mbar, ToF analyser pressure 9.01 e" 7 mbar, 
backing pressure 4.0 e° mbar. All mass spectra in this study were externally calibrated 
using a solution of cesium iodide (10 mg/ml in water) and were processed using 
Masslynx 4.1 software (Waters, UK). 

Dissociation constants. Dissociation constants were derived from the mass 
spectrometry data using the method of Rose et al, (201 1) 27 . Briefly, the intensity, as 
determined by the area under the peak, of all peaks relating to SOD1 monomers (MJ 
and dimers (DJ were summed to give total signal intensity (Its) : 

Im + Id = Its (1) 

The proportion of monomer signal (P M ) was determined by dividing the intensity of 
the monomeric species (I M ) by the total signal intensity (Its): 

Im/Its = Pm (2) 

The concentration of monomer at equilibrium ( [M] eq ) was determined by 
multiplying the proportion of monomer signal (P M ) by the total protein 
concentration in uM as determined by UV/Vis spectroscopy ([P Q ]): 



P M x[P 0 ] = [M] eq (3) 
The concentration of dimer at equilibrium ([D] eq ) is defined as: 

([Po]-N eq )/2 = [D] eq (4) 

The K D was derived from the gradient of the plot of [D] eq against [M] 2 eq (fitting 
performed with Prism 5.0, GraphPad Software, Inc.) In each case at least 5 data points 
were plotted with repeat experiments being performed if the line of best fit gave least- 
squares-regression, R 2 < 0.75. 

Statistical analysis. Statistical analyses were carried out using Prism 5.0 software 
(Graphpad). All mass analyses were performed in triplicate (n = 3) with error bars 
representing the standard error of the mean. The dissociation plots (Figure 2) were 
fitted with a Log (Gaussian) function. The K D plots (Fig. S4) were fitted using linear 
regression. 
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